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Abstract 

00 ' 

We study the risk premium impact in the Perturbative Black Scholes model. The Perturbative Black Scholes 
^ | model, developed by Scotti, is a subjective volatility model based on the classical Black Scholes one, where the 

^s^j . volatility used by the trader is an estimation of the market one and contains measurement errors. In this article 

we analyze the correction to the pricing formulas due to the presence of an underlying drift different from the 
risk free return. We prove that, under some hypothesis on the parameters, if the asset price is a sub-martingale 
under historical probability, then the implied volatility presents a skewed structure, and the position of the 
minimum depends on the risk premium A. 
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1 Introduction 



It is now common knowledge that the Black and Scholes model, which worked well before the 1987 crash, is 
nowadays unable to price options correctly. As can be deduced by comparing the two papers by Rubinstein |23| 
and Jackwerth and Rubinstein [18] . something has evidently changed in the market for options after that event. 
The most shared explanation for the failure in which BS model incurs today is usually thought to reside in the 
\ fact that the constant volatility parameter it proposes is not a good representation of reality anymore. Empirical 
7-H ■ evidence shows that the underlying stock volatility, for example, is not time invariant during the life of an option. 

Moreover, while before 1987 the lognormal distribution of stock prices implied by the Black and Scholes model 

■ seemed to be a good approximation of the real one and volatility observed across strike prices had a moderately 
| pronounced smile, from that date onwards the implied volatility curve appears to be steeper and generally skewed 

■ to the left. Jackwerth and Rubinstein [18] , recovering stock price distribution from observed prices, empirically find 
| a "fatter" left tail phenomenon. Constantinides et al. [8], in the context of an equilibrium model, find stochastic 

■ dominance violations on both tails of the implied volatility curve. Christensens and Prabahla [7], for example, 
00 ] suggest that a regime switch has occurred after the crash. 

The most natural solution to the problem of pricing options more correctly seems then to let volatility change 
with time. Many pricing models have then been proposed, with different formulations for the stochastic process 
driving volatility, e.g. local volatility models, see Dupire [TT], or stochastic volatility models. Then, an option 
pricing model is characterized by a system of differential equations, since two different processes are specified, one 
jf^ ■ for the stock price and one for its underlying volatility. 

The first authors to solve the problem of pricing options with stochastic volatility were Hull and White [T7] . 
However, they were able to obtain closed form solutions only for the case of uncorrelated volatilities and stock 
prices, while Heston [16] . using a new technique, managed to find exact prices also for the correlated case. In this 
stream of literature, one of the most used model in practice is probably the SABR one, introduced in Hagan et al. 
|14j . which provides excellent fitting for interest rates derivatives. More sophisticated models include for example 
the possibility of jumps in the stock price evolution, e.g. see Brigo et al. [5], or directly in the process for volatility, 
see Eraker et al. [10] , 
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However, stochastic volatility models like the ones described above are usually complex and characterized by 
a large number of parameters and, unless in special cases (SABR model, Heston [16]), they do not provide closed 
form solutions for vanilla options prices. 

The PBS model (Scotti [H]) introduced a new category of stochastic volatility pricing models, being founded 
on the notion of subjective volatility. Using error theory through Dirichlet forms, the PBS model generalizes the 
standard Black and Scholes one, imposing an error structure on volatility. Thus, rather than specifying a possible 
pattern of evolution for volatility through time, the perturbative approach deals with the concept of measurement 
errors present in the estimation procedure performed by the trader. 

One of the most important issues of this model lies in the possibility of obtaining closed forms solutions for 
European vanilla option prices and for each kind of derivative which has a closed form solution using the classical 
Black and Scholes model such as Asian and barrier options. This important framework, joint with the flexibility 
of the model, permits us to calibrate it to different markets and fit them, even if they imply opposite behaviours 
of the implied volatility curve. PBS can reproduce a right-tailed or a left-tailed skewness effect, as well as sharper 
or flatter slopes, obviously depending on the calibration of parameters. In Scotti |24j . sufficient conditions for the 
presence of a smile are derived in the case with no drift term in the stock price dynamics. 

This article studies a natural implication of the PBS model: the dependence of option prices on risk premia. 
Dependence of option prices on the expected excess return on the stock is ruled out in the classical Black and 
Scholes model, as a result of the lognormality assumption. Lo and Wang [20], starting from the evidence that 
the predictability naturally captured by the expected return on stocks is affects option prices, assume an Orstein- 
Uhlcnback process for stock prices and show that the impact of a drift term is not negligible anymore. 

In the PBS model we maintain the assumption of lognormality of stock prices. However, we show that the 
presence of measurement errors in the estimation of parameters induces market incompleteness and lets the hedging 
position of a trader not invariant to the expected excess return on stock. We analyze the implications of this fact 
and we show that taking into account the impact of a risk premium on stocks it is possible to reproduce the most 
commonly observed behaviours of the market for options in terms of implied volatility. The PBS model is able to 
generate the usual smile and skew effects pointed out by the empirical literature we addressed above. 

Section 2 briefly summarizes the basic concepts of error theory based on Dirichlet forms, section 3 recalls the 
most important features of the PBS model, section 4 studies the impact of the drift term on the general profit and 
loss function and in the particular case of the price of a European call option, section 5 analyzes the sensitivity 
of the volatility implied by the model to some parameters and above all to risk premium. Finally, the appendix 
provides technical computations. 



2 Preliminaries 

This section resumes the notations used throughout the paper and briefly surveys the key notions of error theory 
using Dirichlet forms0 

2.1 Notation 

We use the following notation: 

• (CI, T ', P) is the historical probability space, sometimes simply denoted with Q. 

• {•^ ? t}o<t<T 1S a filtration defined on the probability space. 

• {^t}o<t<T ^ s the associated Brownian motion, i.e. a Brownian motion adapted to the filtration {^Ft}o<t<T- 



• I fl, T , P) is another probability space, used to represent the uncertainty on the volatility parameter, denoted 



simply with fl. 

• E[»] and E[»|J r t ] denote, respectively, the expectation and the conditional expectation under the probability 
measure P, while E[»] denotes the expectation under P. 

For a complete and exhaustive treatment of error theory and Dirichlet forms, see Bouleau [2], [I], [2] and Fukushima |13) . 
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• (-Pt) t>0 denotes a strongly continuous contraction semi-group. A its generator, with domain T>A, and T the 
"carre du champ" operator associated with the Dirichlct form of the semi-group, with domain D. 

• T denotes the maturity, a fixed positive number. 

• we suppose that there are two traded assets in the market: a stock, with price St at time t, and a risk free 
money account. For the sake of simplicity we assume that the risk free interest rate is zero. 

• (To is the stock market volatility, supposed to be constant, while a is the volatility estimated by the trader 
which we assume to be stochastic. 

• [i is the the expected rate of return on the stock, A represents the risk premium and C the cumulated risk 
premium, i.e. c = \Vt 

2.2 Error Theory using Dirichlet Forms 

Asset pricing models are characterized by a set of parameters that permits to reproduce the real world features 
with some degree of freedom. These parameters have to be estimated in order to be used in practice. The classical 
way to find these estimations is to use market data to make the model fit the prices observed on the market. Hence, 
the final result of this procedure is that we obtain a mean value for each parameter along with some uncertainty 
due to estimation errors. 

What is the impact of such an uncertainty? Generally, models' users consider the mean values of the estimations 
of their parameters and implicitly consider them as deterministic, since their variance is generally small with 
respect to their mean. However, when non linear functions are then used in pricing procedures, the impact of this 
uncertainty, even if small, could be important, since the moments of random variables are distorted. What is the 
amplitude of the bias induced by forgetting the probabilistic nature of an estimated parameter? 

Giving an answer to this question is not an easy task. If we consider a non linear function F and a random 
variable a it is well know by Jensen's inequality that the expected value of F(a) is different from the function 
F evaluated at the expected value of a. A straightforward computation of this difference, however, is usually 
impossible, since the functional forms involved are often difficult to treat. This is the reason why the probabilistic 
uncertainty on the estimation of parameters, whose presence was already known by Gauss, is frequently neglected. 

If a direct approach to solve this problem is unsuccessful, we can try to find a way to reach a solution addressing 
the question from another perspective. We can use the fact that the variance of an estimated parameter is very 
small compared to its mean in order to justify a Taylor expansion; studying the bias and the variance of F(a) we 



where Co is the estimated value of the random variable a, e Bias [a] is its bias, i.e. the difference between <jq and 
E[er], and e Variance [a] is its variance. 

Remark 2.1 The first relation in \2.1\) summarizes the key difference between a deterministic uncertainty and a 
probabilistic one, since the mean of a function evaluated on a random variable is distorted if it is non linear. 

The starting point of error theory using Dirichlct forms consists in considering a very small e and to stop the 
Taylor expansion at the first order. The theory involves thus the use of two operators, a bias operator A and a 
variance-covariance one T, linked together by the chain rule in equation (|2.1j) . Such an environment can be studied 
through semi groups theory. The operator A is the generator of the semi-group and V is the " carre du champ" 
associated with the Dirichlet form of the semi-group. 

The most important framework of error theory is the error structure. We recall its definition: 

Definition 2.1 (Error structure) An error structure is a term 



find: 



(2.1) 



E [F(a) - F(a ) 
E \(F(a) - F(o- Q )) 2 



e {F'(o- ) Bias [a] + \ F"(a ) Variance [a]} + o(e) 
e (F'(o-q)) Variance [a] + o(e) 




where 
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1. (il, J 7 , P^ is a probability space; 

2. D is a dense sub-vector space of L 2 (q, J-, , 



3. r is a positive symmetric bilinear application from D x D into L 1 PJ satisfying the functional calculus 

requirements on the class C x C\Lip. This means that, if F and G are two functions of class C 1 and Lipschitzian, 
u and v £ D, then F(u) and G(v) £ D and 

T [F(u), G(v)] = F'(u)G'(v)T[u, v] P a.s.; 

4- the bilinear form £[u, v] = — E [T[u, v]] is closed; 

5. The constant function 1 belongs to U>, i. e. the error structure is Markovian. 

Hypotheses 2, 3 and 4 together ensure that £ is a Dirichlet form, with T as carre du champ operator. 

We use the simplified notation T[u] = T[u, u] to indicate that the operator T is applied twice on the same 
argument. The couple (r, P) defines a unique semi group (Pt)t>o an d its generator A thanks to the Hille-Yosida 
theorem E|. Hypothesis 5 provides that the semi-group (Pt) t >o is Markovian. 

Therefore, we defined two operators T and A that satisfy the chain rule ()2.1|) . 
It is now useful to conclude this section providing an example of an error structure: 

Example 2.1 (Orstein-Uhlenbeck structure) 

(fi, T, P, D, r) = (k, B(R), ix,H r {ij),T[u, u] = {u'f 



where B(JBL) is the Borel a -field ofM., [i is a gaussian measure and H 1 ^) is the first Sobolev space with respect to 
the measure [i, i.e. u G H l {n) if ' u € L 2 {n) and u' belongs to L 2 (fi) in distribution sense. 
The associated generator has the following domain: 

: u" — x f belongs to L 2 (//)} in distribution sense 

and the generator operator is 

A[u] = \u"-\l-u> 

where I is the identity map on R. 

This example gives the basic idea of an error structure on a parameter. Moreover, as shown in Hamza |15j . it 
is important since every Dirichlet form on R has a characterization. 



3 Perturbative Black Scholes model 

In this section, we recall the key features of the Perturbative Black Scholes model introduced by Scotti [2~i] . 

The PBS model is based on the classical Black Scholes model, (Black and Scholes pQ). If we assume that the 
interest rate is worth zero or, from an economic point of view, that all assets are priced in terms of the money 
market, then the underlying stock price follows the SDE 

(3.1) dS t = [iStdt + ao S t dW t 

where pi is the return on the stock, o~q is the volatility and Wt is a Brownian motion. 

2 see Fukushima 1131 for a complete proof 
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In the BS model pricing formulas depend on the diffusion term only and not on ^; we find closed forms 
expressions for the premium and the greeks of vanilla options. In contrast with its simplicity, unluckily the BS 
model cannot reproduce the so called smile effect: the volatility implied by the BS model is constant across strike 
prices, while the observed one is usually u-shaped. 

We introduce the notation of risk premium A as the ratio between the expected return on the stock and market 
volatility. 



The PBS model lies on three main hypotheses: 

1. the stock price follows a geometrical brownian motion with fixed and non perturbed volatility cto; 

2. the trader has to estimate the volatility parameter, and the value of his own estimation contains intrinsic 
inaccuracies. The model reproduces this fact through an error structure; nonetheless we assume that the 
stock price St is not erroneous. We evaluate the impact of the perturbation generated by those measurement 
errors on the profit and loss process used by trader to hedge a position on a vanilla option; 

3. the trader knows the existence of the perturbation described above and wants to modify his own offered 
prices in order to take into account the bias present on volatility and, as a consequence, on the hedged 



Summarizing, all traders use a geometric Brownian motion to model the stock price process and they hold 
some positions involving vanilla options; they use observed market prices to determine the values of parameters by 
inversion of pricing formulas. Thus, they find an observed volatility process usually known as implied volatility, 
they take it as a forecast for future volatility and hedge their portfolio accordingly. 

Since we made as an assumption that the trader knows the existence of errors in his estimation procedure, 
volatility is incorporated into the model in two different ways. It has a "market" value, the classical parameter 
used to set up standard pricing formulas, and a subjective one. The former is denoted with uo and it is supposed 
to be a constant parameter as in the Black Scholes model. The "subjective" volatility notion comes from the 
intuition that in the real world, when an operator deals with the problem of option pricing, he docs not know 
the precise value volatility will assume during its life. Hence, he has to estimate it from market observations. 
The value he gets from this procedure, as pointed out above, will obviously be subject to measurement erroriH, 
captured in the PBS model by the error structure form. We assume that the stock volatility "market" value is 
also the mean value of volatility in the erroneous estimation procedure performed by the trader. The volatility 
estimated by the trader is then a random variable, and is "subjective", since it can assume a different value in the 
expectation of each operator. 

The profit and loss process of a trader has a key role in the PBS model. The value of this process at maturity 
is given by: 



where F(^o, So, 0) is the security premium, the integral term represents the hedging strategy, $(Sr) is the 
Payoff and S t follows Black Scholes SDE (|3~T|) . 

4 Impact of the drift term in security pricing 

In this section, we study the impact of a non zero drift term in the diffusive process assumed for stock prices on 
prices determined with the PBS model. 

As we have shown previously, the expected profit and loss function from the hedging position is then in turn 
a random variable, characterized by a bias and a variance term, which make it different from the one implied by 
the BS model. We make an important remark: 

3 Measurement errors arise from the uncertainty expected using the central limit theorem. 




position. 



(3.2) 
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Remark 4.1 (Drift impact) In the PBS model, the profit and loss process defined in equation iS. 2\) depends 
crucially on the drift rate fi, which is the expected excess return on the stock. As a matter of fact, the integral term 
in i3.2\) depends on the diffusive process described in \3.1\) where \i plays a role. 

In the Black Scholes model, instead, the price of an option does not depend on the drift term. In that case, in 
fact, the PkL process is worth zero almost surely; as a consequence, we can change the probability measure without 
altering the result. If, as in the PBS model, we assume that the volatility ft used by the trader is a and not o~$, the 
profit and loss process is not worth zero a.s.; on the contrary, it becomes a stochastic process characterized by two 
random sources: 

• the Brownian motion which describes the evolution of the stock price and 

• the process ft, the trader's volatility, which depends on an independent probability space. 

As a consequence, we cannot change the probability measure without changing the value of the profit and loss 
process at maturity. 

We suppose that the trader's volatility ft is the time independent random variable a we defined in the previous 
section. Using the language of Dirichlct forms, we derive the following expansion for the volatility estimation: 

eA[a](a ) + ^eT[a}{a )Af 

where Af is a standard Gaussian random variable. Moreover, we assume that this error structure admits a sharp 
operator. 

We estimate the variance and bias of the error on E [PSzL] . In the computation we assume that a = cto is the 
right value of the random variable, in the sense that if ft = o~o, then P&L(cto) = almost surely. Notice that, 
however, this does not mean that the trader believes the BS model to be correct. 

Then we can prove, see Scotti [24j . that we have the following bias and variance terms: 



AmPkL}} 



(4.1) 




(ctq, S s , s) dS s 



A[<t](<t ) 



da 2 dx 



(ct , S s , s) dS s 



T[ct](ct ) 



T [E [PkL]] = 



■ E 



d 2 F 

da dx 



-I \ 2 



(ct , S s , s)dS s 



r[CT](CT ) 



These values represent the inaccuracies that the trader knows to be present in his estimates. 
We can give the following interpretation to this error structure: 

• the bias in the PSzL process represents a deviation in security prices asked by the trader to the buyer. 

• the variance of the PkL process naturally generates a bid/ask spread on security prices. The width of the 
bid-ask spread depends both on the traders' risk aversion and on the perceived uncertainty on volatility. 

As a consequence of the presence of the error structure, the price of a security is thus not uniquc@, but it can 
be represented, at each instant in time, as a distribution, whose characteristics depend on the parameters which 
characterize the error structure^. 

Therefore, we have shown that the trader must modify his prices in order to take into account the two previous 
effects, namely the variance and the bias on his expected profit and loss process. Thus, he fixes a supportable risk 
probability a < 0.5 and accepts to buy the option at a certain price 



4 When, as in the classical B-S formulation, prices are unique, risk-neutral arguments can be formulated in order to solve the partial 
differential equations which rule the pricing of assets. In this sense, at each instant in time, price can be represented through a Dirac 
distribution. 

5 As pointed out in Scotti (2007), PBS model induces market incompleteness. 
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(Bid Premium) = (BS Premium) + e A [E[P&L]} + ^eT [E[PkL]] N a 

where Af a is the a-quantile of the reduced normal law. Analogously, the trader accepts to sell the option at 
the price 



(Ask Premium) = (BS Premium) + e A [E[PkL]} + T [E[PkL]} M-a 
Since M a + A/i- Q = 0; the mid-premium is 

(4.2) (Mid Premium) = (BS Premium) + e A [E[PkL\] 

and the bid-ask spread is 



(4.3) Bid- Ask spread =2 v / eT [E[P&L]} N a 

4.1 European Call options 

We now focus our attention on European call options and we study the bias and its derivatives in order to derive 
some sufficient conditions for the presence of a smiled behaviour on implied volatility. We know the premium of a 
call option (see Lamberton et al. |19j ) with strike K, spot price x, volatility cto and maturity T, and we know its 
hedging strategy in the usual Black Scholes setting: 



F(a , x, 0) = xM{dx) - KM{d 2 ) 
dF 

Delta = — (a , x, 0) =Af(d 1 ) 
ox 

lnx-\nK + ^-T r - 

where d\ = -= — - — and d 2 = d\ — oqvT. 

<jqVT 



The following results arc classical (see [19]): 



dF , f=e~^ 

-((To, X, 0) = XV 1 - 



d(J() ' V27T 



4.4 — =- ((To, X, 0) = A^l rf 2 

c?erg c V27T 

d 2 F , , x e~5 d i 



r(>o, x, 0) 



and we can easily prove that 
d 2 F , 



(4.5) 



da dx V2tt <j 

d 3 F d^Ss, s)+d 2 (S 8 , S )-d 1 (S a ,s)4(S a , s) _x d?(Ss)S) 
— -y- -(ct , b a , s) = e » iv ; 



where rfi(5's, s) 



InS. - In if + ^(T-s) 

(Tq y/T - S 



and d2(S s , s) = di(S s , s) — ao vT — s 

We apply the Perturbative Black Scholes model to find the corrections it imposes on the expected profit and 
loss process for a trader who is hedging a short position on a plain vanilla European call option. Then the bias on 
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the call premium is given by two terms. The first one is the bias when (i = 0. This case is accurately studied in 
[241. 



(4.6) 



Afj, = Q [C] |cr=ero — X 



e 2 a i 
V2^ 



A 



di d 2 



2a VT 



Now, if we assume that the drift term of the stock price process is non zero, this correction is not sufficient 
in order to hedge the position correctly. We have to study another term, which is the correction when fi ^ 0. 
While when [i = the stochastic integrals in equation (|4.1| are martingales, if fi > we have to evaluate their 
expectations: 



A 



correction 



[C]\„=c 



= E 



T f)2 



d 2 F 



-E 



o da dx 



(er , S s , s) dS s 



A[a]{a ) 



o da 2 dx 



[a , S s , s) dS s 



In appendix [^1 we compute the two integrals and we find 



-4, 



correction 



'C]L = 



■ N{d 2+ C)-Njd 2 ) 1 i4 



K 



(4.7) 



<7 VT 

1 



c 

a VT 



'lit 



aVT 



2tt 
1 



d z 2 + 4 - 2 + — + a VT [ d 2 - - - 



2 7T 



c c 



op VT ( da . 
C { £ C 2 



8 R d 2 



4 d 2 



[Af(d 2 + C)-Af(d 2 )] 



C C 2 



-i d 2 

e 22 



where £ is the cumulated risk premium: 



(4.8) 



C = XVT = 



a VT 



Remark 4.2 We remark that the first correction derives from the bias of the option price. It is then an 

uncertainty coming from the error on estimating the value of volatility. 

The second correction, |^.7[ ) is a consequence of the presence of a bias on the strategy, which introduces uncer- 
tainty on the hedging procedure also. 



It is easy to compute the value of the variance term of the error structure for a call option: 



(4.9) 



T[Call] 



" -.e-^+K 



'2-k 



M(d 2 + C) -N(d 2 ) 
C 



'2 it 



-g 2 "2 



T[aVT] 



We assume that the ask price is then simply the mid price increased by a standard deviation, symmetrically 
the bid price is the mid price decreased by a standard deviation. The spread is simply given by 



(4.10) 



Bid- Ask spread [Call] = 2 Ve T[Call] 
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5 Numerical Analysis 



In this section we explore the sensitivity of the PBS model to some parameters and, through numerical analysis, 
we give evidence of the fact that the model is able to reproduce all the observed behaviours of the implied volatility 
curve. The existence of closed form solutions to the pricing formulas allows us to make some comparative static 
exercises in order to analyze the dependence of the implied volatility curve on time horizon and, above all, on the 
drift term. 

5.1 Parameter sensitivity 

First of all, we point out that the corrections obtained with respect to Black and Scholes prices depend on the 
choice of the parameters A and T and on the magnitude of the e tern0. A captures the bias introduced on the 
profit and loss function, while T is a variance term, e is just a scale factor. It must be small enough to make the 
higher order expansion terms be negligible. Let us analyze the sensitivity of the volatility implied by the PBS 
model to the choice of these parameters. Thus, we fix a value for a and the other real world parameters and we 
let the coefficients of the error structure vary. 

Figure [T] shows the effect of an increase in the absolute value of the coefficient of the bias term. As we will 
explain below, there are good reasons for considering a negative bias. Then, the lower the coefficient, the more 
the curve shifts downwards and the point of minimum variance to the right. Increasing the value of the coefficient 
of the variance term, instead, clearly "opens" up the smile, which becomes more pronounced. Moreover, as can 
be seen in figure a higher coefficient is associated with a more pronounced skew effect which makes the implied 
volatility higher for out-of-the money options, compared to in-the-money ones. 

A change in cpsilon, instead, combines the two effects described above. Notice that this parameter must be 
small enough to justify expansion (|2.f \ in our case. A higher epsilon, thus, produces both a downward shift of the 
curve and more pronounced smile and skew effects, see figure [31 

From now on, we fix the values of these three parameters to make some comparative static analysis of the 
parameters that capture the real world features. We set epsilon to 0.02 for convenience. 

The coefficient on the variance term is set, by a normalization argument^, to 

(5.1) Y[o\ aa =oJ 
implying that 

(5.2) Variance = ea a 2 
The bias coefficient is instead set to 

(5.3) A[a]\ ao =-5a 
leading to 

(5.4) Bias = — 5s<tq 

This last choice is made in order to reproduce a precautionary effect. The hypothesis we make is that for 
some reason, the trader believes he overestimated volatility in his procedure. This feeling can be justified by two 
reasons, one mathematical and one economic. 

6 Since we have three parameters, a possible way to calibrate the model to the market behaviour is by using instruments which 
price variance (e.g. variances swaps). Such derivatives permit to find an implicit link between the bias and the variance term; by 
fixing an epsilon, it is possible then to calibrate the model on just one parameter, see Scotti 1251 . Another possibility is instead to fix 
an arbitrary epsilon small and use the implied spread to calibrate the two coefficients A and T. 

7 Obviously, this coefficient can not be negative. As shown in Figure [2] higher values of this coefficient lead to more pronounced 
smiles. 
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The mathematical explanation lies in the analysis of the usual formula used to estimate historical volatility 
under the hypothesis of lognormality of stock prices. Since it is a concave function, it is more likely that the 
approximated value found by the estimation procedure is an overestimation of the true one. 

The economic explanation lies in the way volatility is usually described in models. Markets are opened for 8 
hours a day only. However, the flow of information does not stop when markets are closed: variability accumulates 
even if securities are not traded. Then, in almost every pricing model, volatility is described as a continuous 
process. This is of course a simplification, but seems nevertheless reasonable. However, it has been shown by some 
authors, see Stoll and Whaley [55] that overnight volatility is consistently lower than intra day one. Hence, it is 
straightforward to believe that usual models overestimate volatility. 

Let us consider the PBS model prediction on a one-month European call option with the parameters we set 
above. 

First, we keep the risk premium measure fixed and we analyze how implied volatility changes with different 
maturities. Our finding is that we obtain curves which are flatter as long as the option time horizon becomes 
longer. This behaviour is consistent with empirical evidence on almost every derivative market, see Hagan et al. 
|14j . Figure 4 shows that the implied volatility curve is skewed to the left for each option; the point of minimum 
variability shifts towards higher strikes for longer time horizons. 

As we showed in the previous section, the use of our perturbative approach implies that option prices and, 
thus, implied volatilities are affected by changes in the risk premium. If we let A , the risk premium we defined 
previously, change and we fix the parameters that characterize the error structure, we can observe and analyze 
this sensitivity. 

Figures and [5] show the behavior of the implied volatility curve on a 1 month European call option for an 
expected excess return on stock term that ranges from to 0.2. The curve evidently shifts to the right side of the 
graph as the risk premium term increases. With almost every value up to 0.2, i.e A = 1, there is a skew effect 
towards lower maturities. The lower the risk premium, the higher is the value of implied volatility for deep in the 
money options and the lower for options which are far out of the money. As the value of [i increases the curve 
appears to become steeper on the right side. In particular, for this parameter choice, for a very high risk premium, 
there is a slight tendency to change the skew directior|§ Curves cross approximately at the money, around 102. 

These findings are consistent with those obtained by other stochastic volatility models, such as the SABR one 
(Hagan et al. [2]). The authors of that model, fitting it on prices of Eurodollar options, obtained those behaviours 
of the implied volatility curve, under the hypothesis that asset prices and volatilities are correlated. 

5.2 Spread Analysis 

Up to now, we have considered the mid price only. As pointed out before, the perturbative approach used by the 
PBS model can naturally generate a spread on prices and volatilities. 

The spread on implied volatility is then obtained by inversion of the pricing formula. 

For the same set of parameters described above, we can thus analyze the effect of changing the drift term on 
a theoretical bid-ask spread. Figures [7] and [5] give an example of price and volatility spread behavior for a chosen 
value of jj,. It is straightforward to notice that higher prices imply higher volatility. Figure [H] shows the magnitude 
of the spread on implied volatility for three different values of fx. For low strikes, the spread is higher when there 
is no risk premium; it reaches a minimum around the money, then it starts increasing. For out of the money 
options, the behavior is reversed: the spread is higher the higher the risk premium. The main difference we find 
with the standard fj, = setting is that the spread has no longer its point of minimum variance around the money. 
The variance spread becomes indeed wider as the strike increases. As shown in figure [TIB the relative spread on 
prices (spread-mid price) is almost zero for in deep in the money options, then increases sharply with both strike 
and risk premium for out of the money calls. 

Let us finally consider directly the effect of the addition of the correction in equation ()4.7j) to the PBS model 
implied volatility curve. Figure [TJJ clearly shows that the presence of a risk premium skews the implied volatility 
curve toward higher strikes. 



8 For fi = 0.2 the implied volatility at moneyncss 1.15 is slightly higher than at 0.85. Unreported simulations show that this 
behavior is common to every choice of parameter. This could suggest that in periods of high risk premia, volatility should tend to be 
higher for out of the money options. 

9 Notice that the analysis of the relative spread leads to the same conclusions. 
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Figure 4: Skewed implied volatility depending on maturity. 
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Figure 5: Implied volatility depending on A. 
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Figure 6: Implied volatility depending on A. 
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Figure 9: Volatility spread depending on A. 
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Figure 10: Implied volatility in PBS model with and without drift. 
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6 Conclusion 

In this article we have studied the Perturbative Black Scholes model, introduced in [21], when we drop the 
hypothesis that the underlying is a martingale under the historical probability. Without imposing any behaviour 
of volatility through time, we showed that the hedging procedure of a trader who estimates it depends on the 
expected excess return on stocks. 

We then introduced a correction with respect to Scotti [H] when the drift term of the diffusive process for the 
stock price is different from the risk free rate. We found a closed form solution for the pricing of a European vanilla 
call option. This formula depends on the same parameters of the classical Black Scholes model, i.e. the volatility 
(To and the parameter d,2, on the two parameters of the PBS model, the variance T[a] and the bias A[a], which 
characterize the error structure of the volatility estimated by the trader. Since the PBS model induces market 
incompleteness, pricing formulas depend also on the cumulated risk premium C, as shown by equation (|4.7|) . 

We analyzed how a simple risk aversion argument forces the underlying price to be a sub-martingale and we 
studied the dependence of implied volatility on the parameters of the model. We numerically studied the most 
case in which the volatility used by traders is an overestimation of the true value and we showed that higher risk 
prcmia tend to increase the skewness and the smile of the implied volatility curve, since the distribution of stock 
prices at maturity is shifted towards higher values. 

We finally found out that the Perturbative Black Scholes model with drift can reproduce the behaviour of the 
implied volatility curve after the 1987 crash. 
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A Computation 



We have to compute 
(A.l) 



t pa. 



8 2 F 
da dx 



(ct , S s , s) dS s 



where St follows the Black Scholes diffusion (|3.1j) and F((Tq, S s , s) is the price of a call option with strike K, 
starting at time s, when the spot value is S s and the volatility is ctq. 
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We integrate by part and we find 
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where A" is the cumulated distribution function of a reduced gaussian random variable. 
The second term that we have to compute in equation (|4.7|) is 
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We can compute this term following the same steps we used for the first term (|A.1|) . 
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Finally, we integrate by parts three times and we find 
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